Experimental evidence for nematic order of cuprates in relation to lattice structure 
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Experiments that have been interpreted as providing evidence that the pseudogap phase in 
cuprates is an electronic nematic are discussed hom the point of view of lattice structure. We 
conclude that existing experiments are not sufficient to prove a nematic origin of the pseudogap in 
underdoped cuprates. 



PACS numbers: 74.72. Kf 



I. INTRODUCTION 



Understanding the nature of the pseudogap^"— and its 
relationship to superconductivity has emerged as a cen- 
tral problem in high Tc cuprate superconductivity. While 
there has been considerable progress in elucidating its 
momentum and doping dependence, its microscopic ori- 
gin and relationship to superconductivity remain contro- 
versial. Some widely discussed experimentally supported 
possibilities are (1) that it is associated with an order 
with a quantum critical point below the superconducting 
dome that may play a key role in superconductivity and 
normal state transporti^ii (2) that it represents Cooper 
pairing without phase coherence, i.e. a strong coupling 
precursor of superconductivity,-^ and in contrast (3) that 
it is a phase associated with cuprates that competes with 
superconductivityi^ii^ There have also been recent dis- 
cussions, also based on experimental results, of the pseu- 
dogap phase as an electronic nematic, specifically asso- 
ciating the pseudogap temperature T* with an onset of 
electronic nematicityJ^^— The purpose of this paper is 
to discuss the experimental evidence for this. We argue 
that, contrary to recent suggestions fii the evidence for 
an association of T* with an onset of nematic order is 
actually weak. 

We emphasize that we are not discussing the so-called 
stripes in the strongly underdoped 214 cuprates. The- 
oretically, these are a dynamic phase separation that 
minimizes the Coulomb energy, and they occur on a 
larger length scale that the microscopic nematicity that 
is thought to occur in Sr3Ru207 and above Tafm in 
pnictides. 

Specifically, we analyze the statement that at the tem- 
perature T* where a pseudogap in the excitation spec- 
trum develops the symmetry the Cu02 plane(s) breaks 
from a four-fold (tetragonal) to at most a two-fold. The 
experiments that seem to support this idea include trans- 
port measurements on YBa2Cu307_5,^^ neutron scatter- 
ing experiments on YBa2Cu306.43^ and especially scan- 
ning probe measurements on several materialsJ^^. 

We begin with definitions, which are important in this 
area because research on nematicity in strongly corre- 
lated systems is quickly developing field and perhaps as a 



result the term has become broadly applied in some cases 
to phenomena that are rather common and do not fit un- 
der the standard definitions. The formal definition of a 
nematic phase is found in the theory of liquid crystals J^ 
There a mesomorphic (lacking long range order) liquid 
crystal is called nematic if the constituent anisotropic 
molecules have one preferential orientation axis. At the 
same time, the centers of molecules are disordered. In 
this sense, a nematic solid state phase would be a glass or 
a disordered alloy where short-range correlations in one 
direction are different from the other. This is illustrated 
in Fig. [T] The left panel shows an ordered orthorhombic 
crystal; note that the lattice parameters a = b, but the 
crystal symmetry is orthorhombic and a long range order 
is present. Such as state is not a nematic. Note also that 
the two species need not be different ions, they can just as 
well be spins on different sites. The middle panel shows a 
textured crystal where there is a long range order in one 
direction, but not in the other. Strictly speaking, this is 
not a nematic phase either. On the other hand, the right 
panel is an example of a canonical nematic phase: there 
is no long range order at all, yet the number of the like 
bonds along a is 70% larger than the number of the unlike 
bonds, while along b there is no such correlation. Such a 
crystal may show substantial anisotropy in e.g. transport 
properties. In fact, this is the signature of nematicity - 
a breaking of rotational symmetry without translational 
symmetry breaking. 

Nematicity in the electronic subsystem may or may not 
trigger a significant difference between a and b depending 
on the relative energy scales associated with the nematic 
order and the lattice (i.e. elastic constants). For in- 
stance, a considerable anisotropy in transport is observed 
in the high field nematic phase of Sr3Ru2 07 without any 
detectable a/b differentiation. ■^^ On the other hand, there 
is a strong structural signature at the onset of the very 
likely nematic phase that precedes the spin density wave 
transition as as LaFeAsO is cooled. ^^'^^ However, note 
that the spin density wave phase itself has long range 
translational symmetry breaking, and while it lowers the 
rotational (point group) symmetry from tetragonal to or- 
thorhombic, it is not a nematic. 

To summarize, the concept of electronic nematicity can 
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FIG. 1: Schematic depiction of different types of ordering in 
crystalline media: (a) Symmetry-breaking long range order; 
(b) Symmetry-breaking texture (long range order in one di- 
rection only); (c) nematic order. In the last case, there is 
no long-range ordering in any direction but the number of 
the like bonds along x is different from the number of unlike 
bonds, while along y there are the same number of like and 
unlike bonds. 



and should be applied to crystalline solids to describe a 
situation where electronic or magnetic interactions cause 
the rotational symmetry to be spontaneously broken 
without long range positional or magnetic order. In fact, 
only one truly tetragonal cuprate,(Ca,Na)2Cu02Cl2, has 
been studied in this context, but only by scanning tun- 
neling microscope (STM) based spectroscopy, which is 
a surface sensitive technique, and there is evidence that 
the surface may have a lower symmetry due to a recon- 
struction (see below). Therefore the issue is the more 
subtle one of whether there is an electronic instability of 
the Cu02 planes towards a breaking of their four fold 
symmetry that drives the pseudogap. This must be dis- 
tinguished from the possibility that specific probes are 
particularly sensitive to the pre-existing lack of four-fold 
symmetry of the lattice. This is complicated by the fact 
that electronic states around a gap or pseudogap are gen- 
erally much more sensitive to anisotropy than the more 
dispersive bands that would be present prior to the gap- 
ping. 



II. MATERIALS 

When discussing electronic instabilities in complex ma- 
terials, such as cuprates, it is important to keep in mind 
that electrons are charged particles. Therefore they are 
subject to the Coulomb interaction both with other elec- 
trons and with the ions that make up the crystal lattice. 
The effective strength of the interaction between quasi- 
particles at the crystal lattice is affected by screening. 
In fact, it is generally the case that different electronic 
states and different energies have different couplings to 
the lattice, as seen for example in mode and momentum 
dependent electron-phonon couplings. There is no a pri- 
ori reason for assuming that such couplings are negligible, 
especially on low energy scales such as those associated 
with the pseudogap. 

The shadow Fermi surfaces seen~ in angle re- 
solved photoemission spectroscopy (ARPES) measure- 
ments provide an instructive example. A low inten- 
sity replica of the Fermi surface of superconducting 
Bi2Sr2CaCu20g+5 was observed shifted from the main 



Fermi surface by (7r/a,7r/a), which is the antiferromag- 
netic ordering vector of the undoped cuprates. This 
shadow was ascribed to scattering by antiferromagnetic 
fluctuations, which would be a novel electronic effect 
presumably related to strong correlations. This led to 
a large body of follow on experimental and theoreti- 
cal work aimed at elucidating the detailed underlying 
mechanism. However, it was pointed out early on that 
the shadow bands were remarkably sharp for a fluctu- 
ation induced feature and that the lattices of the Bi- 
cuprates have structural distortions that include the cell 
doubling needed to produce this shadow i^i^ More recent 
experiments showed that this prosaic picture explains the 
shadow Fermi surface<^"— Importantly, this distortion 
lowers the symmetry from tetragonal to orthorhombic 
and as evidenced by the observation of the shadow Fermi 
surface couples significantly to the electronic states in- 
volved in superconductivity and the pseudogap. Another 
instructive example is provided by the high sensitivity of 
superconductivity in the (La,Ba,Sr)2Cu04 system to the 
structural phase transitions in particular the switching 
between low temperature orthorhombic (LTO) to the low 
temperature tetragonal (LTT) distorted structures with 
composition.^^ 

YBa2Cu307_A- has two Cu-0 systems, Cu02 planes 
and CuOi_5 linear chains. The apical O atoms of the 
plane Cu are also nearest neighbors of the chain Cu in this 
structure. This leads to hybridization between the bands 
at the Fermi surface that are derived from the plane Cu 
and electronic states associated with the chains, includ- 
ing near 5—Q, the chain Fermi surfacci^ The fact that 
there is coupling between these two electronic subsystems 
is seen from the relatively low anisotropy of the resistiv- 
ity. Pel Pah^ i-e. the relatively high c-axis conductivity as 
compared to other layered cuprates and also the opening 
of a substantial presumably induced superconducting gap 
on the chain derived Fermi surface. As mentioned, the 
CuOi_5 layer consists of Cu-O-Cu chains. The chain for- 
mation is associated with O structural ordering, which 
takes place at high T and leads to an increasing or- 
thorhombicity as the O increasingly order with decreas- 
ing T and annealingj ^^'^^ This chain formation persists 
to low doping (high S) and in fact at 8 =0.5 an ordered 
orthorhombic structure with alternating fully oxygenated 
chains and empty "chains" forms. High resolution x-ray 
scattering studies have shown that there is an interplay 
between this orthorhombicity and superconductivity, as 
reflected in an anomaly in hja at T^^ Importantly, sub- 
stantial orthorhombicity persists with increasing b up to 
the superconductor-antiferromagnet phase boundaryi^ 

(Ca,Na)2Cu02Cl2 is a true tetragonal material as was 
shown by neutron diffraction measurements.'^'* However, 
the only experimental study pointing at nematicity in 
this material is by STM, which is surface sensitive. A 
low energy electron diffraction (LEED) image of the sur- 
face of this material was reported by Kohsaka and co- 
workers,'^^ They do not find evidence for a reconstruc- 
tion in their pattern. However, based on the intensity 



of the Bragg peaks in the LEED pattern, it was not of 
sufficient quahty to detect a surface reconstruction con- 
sisting of displacements of the surface atoms (see Refs. 
ISa andlSTI for reviews of the LEED technique). Moreover, 
a strong shadow band is seen in ARPES, and as in the 
Bi compound discussed above, it is sharp in momentum 
space. This sharpness shows that it comes from a static 
order rather than fluctuations, which would introduce a 
width related to the momentum distribution of the fluc- 
tuations that provide the scattering. This shadow is then 
a strong indication that the surface is reconstructed. In 
fact, the ionic charges of CI, Ca and Na are almost cer- 
tainly -1, +2 and -1-1, respectively. There is then no 
natural neutral (001) cleavage plane in this crystal and 
so a reconstruction is highly likely. More careful investi- 
gation of this surface is suggested to determine the actual 
reconstruction pattern and its symmetry. 

Finally, it is important to note that (1) an orthorhom- 
bic point group does not contain a four-fold axis and 
therefore regardless of sample orientation, measurement 
direction or coordinate system one cannot have a break- 
ing of four-fold symmetry, since there is no four-fold axis 
to start with; and (2) the symmetry of a surface can 
be lower than the symmetry of the bulk, but it cannot 
be higher; in particular a (001) face of an orthorhombic 
crystal will not have four-fold symmetry. 



III. SCANNING TUNNELING MICROSCOPY 
BASED SPECTROSCOPY 



We start with STM, as these experiments have 
been most emphasized in connection with possi- 
ble nematicity. The systems studied include Bi- 



2212 phases and (Ca,Na)2Cu02Cl; 
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Importantly, 



(Ca,Na)2Cu02Cl2 is the only cuprate that has been stud- 
ied in the context of nematicity at T* and is at the same 
time crystallographically tetragonal in bulk. The mea- 
surements are done at low T, well below T*. Spectra 
are taken as a function of tip position and two energy 
scales are identified in the spectra: Aq, which is related 
to superconductivity and Ai, which is associated with 
the pseudogap. Ai(r) is found to vary considerably with 
position r on the sample surface. This indicates a sen- 
sitivity of electronic states Ai at the pseudogap edges 
to the local chemistry or structure variations on the sur- 
face. This is as expected because states at and near a 
gap edge are in general heavier than those associated 
with dispersive states, and therefore more easily local- 
ized. This helps to confirm the association of Ai with 
the pseudogap. Real space images of the spectra with 
the intensity of Z{r,E), which is the ratio of differential 
conductances and opposite bias, defined in Ref. [l^, also 
show variations, and these are correlated with the varia- 
tions in Ai, so that real space plots of /(r) = Z(r, Ai(r)) 
show strong spatial dependence. What is being plotted 
then is a function that has the most inhomogeneity in real 
space presumably because of it has the most sensitivity 



to variations in local structure and chemistry. 

As emphasized in the STM papers, this function shows 
chain-like correlations running along the directions of the 
Cu-O-Cu bonds with lengths of several atoms. As men- 
tioned the structures of Bi-cuprates have strong lattice 
distortions. They consist of shifts of O atoms in the Bi- 
O planes by ~0.5 A, as well as large mainly transverse 
shifts of the apical O and buckling of the Cu02 planes, 
corresponding to tilts of the CuOe octahedra. Such tilt 
type instabilities come from shifts of the O atoms and 
not the cations. There is also an incommensurate mod- 
ulation associated with the mismatch between the Bi-0 
and Cu02 systemsi^ Further, the Bi^^-O^" surface is 
not charge neutral and therefore additional distortions 
may be expected, although these have not been quanti- 
fied experimentally. In any case, tilt instabilities in per- 
ovskites are well studied4S Layered perovskite zone cor- 
ner lattice instabilities with tilts around [110] are close 
in energy to structures with [100] tilts. Also, because of 
the bonding topology, in particular the square network 
of corner sharing octahedra the tilts of nearby octahedra 
are generally correlated. In other words, clockwise ro- 
tation of an octahedra drives counterclockwise rotation 
of its neighbors because of the shared O atoms. This is 
seen in the phonon dispersions of perovskites, where for 
example in cubic perovskites the R25 branch disperses 
weakly towards the M point but strongly towards F (see 
e.g. Ref. Hll for a discussion of this). 

These tilts, which are strong in Bi-cuprates, provide 
a natural explanation for chain-like correlations seen in 
the STM experiments if (1) local disorder on the surface 
is strong enough to rotate the tilt axes away from [110] 
towards [100] or [010] and (2) this couples to the elec- 
tronic structure being measured. In regard to (1), we 
note that the STM images for the underdoped sample 
show substantial disorder, besides the chain-like correla- 
tions. Visually, this appears stronger in the spectroscopic 
maps that show the electronic structure of the Cu-0 sys- 
tem, and is therefore sensitive to O position than in the 
topographic image (Fig. 4d of Ref. [I^, which shows the 
Bi positions) To quantify this, we took the STM images 
of the spectroscopy Z functions showing the chains (Figs. 
4a, 4b and 4e) of Kohsaka and co-workers^^ and Fourier 
transformed (FT) it. The resulting transforms is shown 
in Fig. [2] Only the first three Bragg reflections (100), 
(110), and (200) are clearly seen and these are broad and 
fall off rapidly with order. Note that even with strong 
variations in the intensities of the STM image from site 
to site, if the atoms were on periodic sites Bragg peaks 
should be seen up to higher order although they may be 
broadened, similar to what is seen in x-ray crystallog- 
raphy for alloys. The FT instead indicates shifts of the 
atoms from the periodically arranged sites in the bulk. 
In regard to (2) we note that the Fermi surface comes 
from Cu dx2_y2 - O Pa anti-bonding states, and the O p 
orbitals {x,y,z) are orthogonal. This is manifested in the 
shape of the Fermi surface which is flat near the nodal di- 
rections. A square Fermi surface can be regarded as the 



intersection of segments of two one dimensional Fermi 
surfaces, and so if one concentrates on the nodal regions 
that provide the gap edge states around Ai, one has the 
intersection of two one dimensional surfaces, one along x 
and the other along y. The a hopping matrix elements 
are reduced quadratically as the bond angle deviates from 
180°, providing a mechanism for coupling of the gap edge 
states to the tilts. As noted, the states at the Fermi sur- 
face are known to couple to the tilts as evidenced by the 
shadow Fermi surface, and focusing on the states at the 
edge of the pseudogap is expected to enhance this sensi- 
tivity as discussed. Therefore, the images shown in the 
STM experiments can have a rather prosaic explanation. 

Returning to chain-like correlations in perovskites, we 
mention that an analogous situation can occur in the lat- 
tice dynamics of the soft mode of ferroelectric perovskites 
such as BaTiOa and KNbOa.— i^ The real space interac- 
tions for the soft mode have a one-dimensional character, 
leading to chain-like correlations seen in diffuse scattering 
experiments above the Curie temperature. These correla- 
tions are of comparable length to the chains in the STM 
images, although unlike them they are dynamic rather 
than static. These are cubic crystals, not nematics, since 
the correlated chains run equally along all three Carte- 
sian directions. Similarly, even though each molecule in 
a liquid of rod-like molecules may have a certain orien- 
tation that correlates with its neighbors it is not a liquid 
crystal unless the global rotational symmetry is lost. 

Examining the FT image (Fig. [2]), one notes that there 
is no clear breaking of x - y symmetry. This means that 
the surface is not obviously nematic on the scale of the 
STM image shown (~30 Cu spacings on a side). One 
may ask if the image is a superposition of small nematic 
domains. In fact, one can find regions in the image where 
two or more neighboring chains appear to run together 
along the same direction. However, one can also find 
regions with a single chain. The FT image shows diffuse 
intensity along the -ffO and QK directions reflective of 
chain-like correlations (see Ref. 14^, but it is not clear 
that other streaks expected are present, nor is there an 
strong asymmetry above the noise level in the shape of 
e.g. the 10 Bragg peak that would indicate local nematic 
domains. We conclude that at least in this image such 
regions, if present above a statistical level, are very small. 
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FIG. 2: (color online) FFT of Figs. 4a (left), 4b (middle) 
and 4e (right), from Ref. ilS . with low order Bragg positions 
indicated. The FFT was done on the figure converted to 
greyscale after trimming a border at the top to remove the 
caption and slivers from the left and right to compensate the 
top border thus yielding square images. 




FIG. 3: (color online) Calculated Fermi surface of 
YBa2Cu306.5 showing the plane derived sections only (see 
text). There are two barrel shaped sections corresponding to 
the two Cu02 planes in the structure and additionally a zone 
folding due to the alternation of filled and empty chains in 
the CuOi_j, 5 = 0.5 layer. 



IV. TRANSPORT 

The Nernst coefficient, i/, is a very sensitive probe of 
gapping.'**''*^ Like the thermopower, it is proportional to 
T at low T and in a simple band structure it is inversely 
proportional to the distance of the band edge from the 
Fermi energy. Unlike the thermopower, v is proportional 
to the mobility, and electron and hole contributions are 
additive. Large Nernst effects are found in semimetals, 
with Bi being a good example,^* while, in normal met- 
als the Nernst coefficient is very small. Therefore, in 
a pseudogapped system one expects the Nernst coeffi- 



cient to grow strongly as regions of the Fermi surface in 
the nodal directions are increasingly removed from par- 
ticipation in transport, consistent with what was found 
experimentally.^"*^ On the level of the standard Boltz- 
mann kinetic theory with a single relaxation rate, the 
Nernst coefficient is defined by derivatives of the trans- 
port function with respect to the chemical potential, sim- 
ilarly to the Seebeck coefficient 4ii^ Thus, the Nernst 
coefficient is particularly sensitive to those parts of the 
Fermi surface that most rapidly change with the energy. 
We show the plane derived sections of the Fermi surface 
of YBa2Cu306.5 in Fig. [31 These were obtained with 



the general potential linearized augmented planewave 
method^ as implemented in the WIEN2k coder^ using 
the generalized gradient approximation of Perdew, Burke 
and Ernzerhof."*^ Our Fermi surface is similar to that re- 
ported previously by Carrington and Yelland.^^ For the 
present discussion, we have removed the chain derived 
Fermi surface from plot and and colored the remaining 
Fermi surfaces according to the local transport function, 
that is, the squared Fermi velocity. We see that despite 
the fact that the Fermi surface geometry retains the ap- 
proximate tetragonal symmetry, the transport function 
appears quite anisotropic, with some area of heavy elec- 
trons at the "y" faces, but not "x" faces. From this plot 
it is clear that hybridization with chain Cu-0 system can 
make transport properties, especially the Nernst coeffi- 
cient, quite anisotropic, even if one excludes the transport 
in the CuO chains themselves. 

Indeed, Ando and co-workers^^ have measured the re- 
sistivity of de-twinned YBa2Cu307_5 crystals as a func- 
tion of oxygen content and temperature. They find that 
the resistivity is anisotropic in all samples, including ^=0, 
and that in underdoped samples ((5=0.55 and (5=0.35) 
the resistivity anisotropy increases with decreasing tem- 
perature. The resistivity for high doping {6=0 and 
(5=0.17) is qualitatively different from the underdoped 
samples, presumably reflecting the fact that the CuOi_5 
chains no longer contribute significantly to the conduc- 
tivity so that, as discussed by Ando and co-workers, the 
anisotropy reflects anisotropy of the Cu02 planes. Im- 
portantly, the anisotropy is a smooth function of tem- 
perature up to the highest reported r=300K, which is 
well above T*, and furthermore shows no clearly appar- 
ent structure that can be associated with the pseudogap, 
although the pseudogap T* can be mapped from analysis 
of the T dependent curvature of the resistivity J^ This is 
consistent with the view that the observed anisotropy is 
a one-electron effect reflecting the anisotropy of the in- 
plane Fermi velocity induced by the hybridization with 
chains. 

On the other hand, the scattering rate may also be 
anisotropic (even if the chain electrons are localized, they 
can still affect the scattering rates). In fact, optical 
conductivity measurements suggest an anisotropy in the 
scattering rate that behaves similarly above and below 
T* 1^ Either way, it is clear that the observed anisotropy 
in v likely arises from anisotropy in the electronic struc- 
ture that exists independent of T* , which manifests itself 
particularly strongly in v, at low T, where the gapping 
is strong. Therefore, while the Nerst anisotropy, may, of 
course, indicate a symmetry breaking, a simpler expla- 
nation in terms of the always present anisotropy of the 
underlying electronic structure, enhanced at low temper- 
ature by gapping is not at all excluded. 

V. NEUTRON SCATTERING 

Inelastic neutron scattering offers a number of advan- 
tages. It is a true bulk sensitive probe and it directly 



measures the dissipative part of the magnetic susceptibil- 
ity, x" &s a function of momentum and energy. Hinkov 
and co-workersi^ reported such measurements for de- 
twinned crystals of underdoped YBa2Cu306.45. They 
showed by analysis of the II,K,L dependent structure 
factors that the magnetic response that they measure 
is from the Cu02 planes and not the chain Cu. As ex- 
pected they find inelastic scattering near the (1/2,1/2) 
position corresponding to the antiferromagnetic ordering 
vector of the undoped compounds. At high energy (50 
meV) the distribution of x" around this momentum is 
isotropic but as the energy is lowered below 7 meV they 
find an anisotropy that increases as the energy is lowered 
towards the quasielastic regime (3 meV). The scatter- 
ing is elongated along the a* direction, indicating the 
presence of chain-like magnetic correlations along the b 
direction. This anisotropy was found to be strongly T 
dependent and vanishes below 100 K. This very inter- 
esting behavior, i.e. the development of ID magnetic 
correlations in the Cu02 planes as T is lowered could 
be an indication of nematicity, or it could be interpreted 
as nearness to magnetism in an anisotropic (orthorhom- 
bic) material. Clearly, further investigation is needed. 
However, one issue that should be emphasized is that 
the anisotropy is already absent when the temperature is 
raised to 100 K, while T* for this composition is at much 
higher T '^250 K.^^ This calls into question the relation- 
ship between this anisotropy and the pseudogap. One 
experiment that would be particularly interesting would 
be inelastic neutron measurements for a truly tetragonal 
material such as (Ca,Na)2Cu02Cl2 when suitable sam- 
ples become available. 



VI. SUMMARY 



In summary, we argue based on symmetry and other 
considerations that the existing experimental data are 
not sufficient to prove a nematic origin of the pseudogap 
in cuprates. Bulk sensitive measurements on tetragonal 
materials, while difficult, will be very helpful in estab- 
lishing whether the phase below T* is associated with 
the onset of nematic order. 
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